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H . Abstract 

An exploratory study is made of the nucleon-nucleon s-wave scattering lengths 
in quenched lattice QCD with the Wilson quark action. The ir-N and tt-tt scattering 
lengths are also calculated for comparison. The calculations are made with heavy 
quarks corresponding to m^/mp ~ 0.73 — 0.95. The results show that the N-N 
system has an attractive force in both spin-singlet and triplet channels, with their 
scattering lengths significantly larger than those for the 7T-./V and tt-tt cases, a trend 
which is qualitatively consistent with the experiment. Problems toward a more 
realistic calculation for light quarks are discussed. 



The understanding of low energy nucleon-nucleon(iV-iV) interactions has been one of 
the most fundamental problems in nuclear physics. A large effort to construct elaborate 
meson exchange models has enabled us to describe the low energy N-N phase shifts and 
deuteron properties quite well0. However, these models are essentially phenomenological 
in that they use mesons and baryons as the basic degrees of freedom and choose the 
coupling strengths among them to fit the experimental data. Deriving the low energy 
properties of the N-N interaction directly from the dynamics of quarks and gluons is 
clearly a problem posed to numerical studies of lattice QCD. 

An important parameter characterizing the low energy N-N interaction is the s-wave 
scattering length. The experimental values are listed in Table 10 together with those 
for 7T-N and tt-tt cases for comparison. We observe that the N-N scattering lengths are 
markedly larger. We also note that the negative sign for the spin triplet channel reflects 
the presence of the deuteron bound state as follows from Levinson's theorem. The small 
values for the ir-N and tt-tt cases are understood as a consequence of chiral symmetry 
for u and d quarks; current algebra and PCAC have successfully predicted the sign and 
magnitude0. Since no such constraint exists for N-N scattering, its large scattering 
lengths are dynamical phenomena of QCD. 

An elegant method for calculating scattering lengths in lattice QCD is provided by 
the formula 0] relating the s-wave scattering length ao between two hadrons h\ and h% 
to the energy shift SE = E^^ — (m^ + mh 2 ) of the two hadron state at zero relative 
momentum confined in a finite spatial box of a size L 3 . It is given by 

6E = -^(1 + c* + cA 2 ) + O(L-t) (1) 

with /i = m hl m h2 /(m hl + m h2 ) and c\ = —2.837297, c 2 = 6.375183. The energy shift and 
hadron masses m^ and m^ can be extracted from hadron four- and two-point functions 
calculated through numerical simulations. Previous studies have shown the practical 
feasibility of the application of the formula for the case of ir-n scattering both for I = 20 
and also for the technically more difficult case of I = 00. With the development of the 



technique to calculate hadron four-point functions^], we have attempted an exploratory 
study of the N-N s-wave scattering lengths, which we report in this article. 

A realistic calculation of the N-N scattering lengths faces three basic obstacles. The 
large values of scattering lengths of order lOfm to be obtained mean that lattice sizes 
much larger than 2 — 3fm, which are accessible in current numerical simulations, will 
be needed to suppress 0(L~ 6 ) corrections neglected in (|1]). A further difficulty in the 
spin triplet channel is that the lowest scattering state orthogonal to the bound deuteron 
state has to be constructed to apply (jl|). Finally statistical fluctuations in N-N four- 
point functions grow rapidly toward large time separations and small quark masses as 
we shall discuss in detail below. In order to avoid these problems in our initial study, 
we have carried out simulations for heavy quarks with m n /m p rs 0.73 — 0.95 in the 
quenched approximation. We employ the Wilson quark action and work at the inverse 
coupling constant (3 = Q/g 2 = 5.7 which corresponds to the lattice spacing a fh 0.14fm. A 
calculation of tt-N and it-it scattering lengths at the same /3 is also carried out to compare 
them with the results for the N-N case. 

The parameters of runs are tabulated in Table 2. Gauge configurations are generated 
for the single plaquette action separated by 1000 sweeps of the pseudo-heat bath algo- 
rithm. Calculations for the it-it and tt-N cases are made on a 12 3 x 20 lattice, while for 
the N-N case we employ a 20 4 lattice anticipating larger scattering lengths. We extract 
the energy shift 5E = E^^ — (m^ + mh 2 ) from the ratio of the hadron Green's functions, 

R(f] <O hl (t)O h2 (t+l)O hl (0)Ok 2 (l)> , , 

U < O hl (t)O hl (0) >< O h2 (t+l)O h2 (l) > { ] 

*=T exp(-6Et) 

where the operators for h% are shifted by one lattice unit in the temporal direction to 
avoid mixing of color-Fierz transformed contributions 0. 

Quark propagators for the N-N case are calculated with a wall source at t = or 1, 
fixing gauge configurations to the Coulomb gauge in all space-time. Projecting out spin 
singlet and triplet combinations of the two nucleon system is made by non-relativistically 



combining the upper components of the nucleon operator iV = { t qC~ l, j5q)q taking account 
of isospin factors. 

Calculation of the tt-N and it-it four-point functions requires quark propagators con- 
necting two arbitrary space-time sites. We handle this problem by the method proposed 
in our previous work0: quark propagators are calculated with a wall source for every 
time slice without gauge fixing. The nucleon source operator for the tt-N case, however, 
is fixed to the Coulomb gauge to increase signal to noise ratio. 

In Fig. 1 we show R(t) for the spin singlet and triplet channels in the A^-A" case 
at K = 0.160 corresponding to m^jm p = 0.85. A clear signal with a positive slope 
is observed for both channels, which means attraction (5E^^ = E^n — 2mN < 0) as 
expected. Similar results are obtained at two other values of the hopping parameter 
K = 0.15(m n /mp = 0.95) and 0.164(m 7r /m p = 0.73). We extract the energy shift SE^n 
by fitting R(t) to a linear form R(t) = Z(l — dE^Nt). The fitting range is chosen to be 
4 < t < 9 for K = 0.150 and 0.160. For the case of K = 0.164 we used 2 < t < 6 for 
the fit due to the poor quality of our data. The fitted values of SE NN are quite small 
(« 0.01), justifying the use of a linear function instead of an exponential. 

For heavy quark the A^-A" interaction becomes shorter ranged since pions exchanged 
between the nucleons are heavy, while the size of the nucleon will not be much reduced. 
In a simple potential model with a linear confining potential, for example, the hadron size 
scales as a third power of the constituent quark mass which is a slowly varying function 
of the hopping parameter. We may expect then that the deuteron becomes unbound for 
heavy quarks. Assuming this to be the case we can extract the scattering lengths through 
([!]) for both spin singlet and triplet channels. The results in lattice units are tabulated in 
Table 3. 

Our results for the n-N and tt-tt four-point functions are plotted in Fig. 2 for K = 
0.164. The data are also fitted with the linear form R(t) = Z(l — 5Et) over 4 < t < 9. 
The resulting scattering lengths are listed in Table 3 (we could not extract the 1=1/2 
tt-N scattering length at K = 0.1665 due to large errors in our data; this point will be 



discussed below). The predictions of current algebra and PCAC are given by 

«o(™) = +-!-*£ (J = 0), -i-*£ (J = 2), (3) 

16tt f£ 8tt /^ 

and 

ao(7rJV) = + ^^" (J= 2 ) ' "^7f (/= 2 } ' (4) 

fi nn = m 7r /2 and /z„.jv = m n m N /(m w + m^r) being the reduced masses. We list in Table 3 
the numerical values of these predictions for the values of m^, m^ and /„- measured in our 
simulations with the tadpole-improved renormalization factor for f n ^\. Our results for 
the scattering lengths are consistent with these values within one standard deviation. It is 
somewhat surprising that the agreement holds even at quite heavy quarks corresponding 
to m^/nip Rj 0.73. ( For the ir-ir scattering lengths a similar agreement has previously 
been reported^], |6|.) 

In Fig. 3 we plot the scattering lengths as a function of m 7T /m p . It is apparent that the 
A^-A" scattering lengths are substantially larger than the tt-N and tt-ti scattering lengths 
already for a heavy quark corresponding to m n /m p ph 0.73. Also noteworthy is the trend, 
albeit with sizable errors, that the values for the spin triplet channel are larger than those 
for the singlet channel, indicating a stronger attraction in the triplet channel, which is 
consistent with the existence of the deuteron bound state. We find these results to be 
encouraging, although the scattering lengths we obtained are still much small compared 
to the experiment: our results correspond to a>o(NN) « 1.0 — 1.5fm in physical units with 
a ~ 0.14fm determined from the p meson mass extrapolated to the chiral limit. These 
small A^-A^ scattering lengths are most likely to arise from a short interaction range for 
heavy pions (m, ~ 0.7 — 1.5GeV) for which our simulations are made; we expect that 
they increase if simulations would be made with a small quark mass|8[]. 

Let us now discuss the issues for extension of the present work toward a more realistic 
case of light quarks. On the theoretical side, if deuteron is a bound state only for light 
enough quarks as assumed here, the scattering length for the triplet channel extracted 
from the lowest two-nucleon energy should diverge at the value of K corresponding to 



the bound state formation, beyond which we need to construct the lowest scattering state 
orthogonal to the deuteron to extract the scattering lengths. Observing the diverging 
trend would be an interesting problem in its own light. 

From the view point of simulations reducing statistical errors presents a major problem 
toward light quarks. A simple argument indicates that the error of R(t) grows exponen- 
tially with t as SR(t) oc exp((2mjv — 3m n )t) for the N-N four-point function; since the 
pion mass vanishes in the chiral limit while the nucleon mass stays non-zero, the slope 
&nn = ^rriN — 3m n becomes larger for lighter quarks. A rapid growth of error is clearly 
seen in Fig. 1; its t dependence is roughly consistent with the estimate. We also expect 
that the magnitude of errors increases for lighter quarks. Thus a progressively larger 
statistics is needed to obtain reliable results as the quark mass is decreased. 

A similar situation holds for the 7T-N case. The estimated slope a n N = ttln — rn n /2 is 
even larger than a^N for the range K = 0.1665 — 0.15 where we have made our simulation 
(e.g., a n N = 0.84 and ajvjv — 0.66 for K = 0.164). This explains a faster increase of errors 
with t for the 7T-N case (compare Fig. 1 and Fig. 2(a)). An increase of magnitude of errors 
is also observed as the quark mass is reduced from K = 0.164 to K — 0.1665. This is the 
reason why we could not obtain the I = 1/2 n-N scattering length for K = 0.1665. 

To summarize, we have explored the possibility of extracting the N-N scattering 
lengths from numerical simulations of lattice QCD. Albeit carried out at heavy quarks, 
our results show that the attractive forces work in both channels of s-wave N-N scattering, 
its scattering lengths being significantly larger compared to the n-N and n-n scattering 
lengths. It remains an important and realistic problem to confirm that reducing the quark 
mass results in an increase of scattering lengths to the values observed experimentally. 
Such a calculation, however, requires a much larger lattice and smaller quark masses than 
are possible at present, and hence has to await a development of the computing power. 

Numerical calculations for the present work have been carried out on HITAC S820/80 
at KEK. This work is supported in part by the Grants-in-Aid of the Ministry of Education 
(Nos. 03640270, 05NP0601, 05640325, 05640363). 
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Tables 



Table 1: s-wave scattering lengths in units of fm| 



N-N 


3 ^1 




-5.432(5) 




*So 




+20.1(4) 


tx-N 


I = 


1/2 


+0.245(4) 




I = 


3/2 


-0.143(6) 


7T-7T 


I = 





+0.37(7) 




I = 


2 


-0.040(17) 



Table 2: Parameters of simulations. All runs are made at /3 = 5.7 in quenched QCD. 



K lattice size # conf. 



N-Ar 


0.150 


20 3 x 20 


20 




0.160 


20 3 x 20 


30 




0.164 


20 3 x 20 


20 


ir-N 


0.164 


12 3 x 20 


60 




0.1665 


12 3 x 20 


30 


71-71 


0.164 


12 3 x 20 


70 



Table 3: Scattering lengths in lattice units at /3 — 5.7 in quenched QCD. Numbers in 
square brackets are current algebra predictions evaluated with the measured values of 
to,,., mjy and /„. with the tadpole-improved Z factor. 



K 




0.150 


0.160 


0.164 


0.1665 


m n 




1.0758(51) 


0.6919(63) 


0.5081(35) 


0.3588(82) 


m p 




1.1302(66) 


0.8171(86) 


0.7008(67) 


0.588(20) 


m N 




1.788(11) 


1.302(13) 


1.093(20) 


0.951(42) 


N-N 


3 ^1 


+ 10.8(1.2) 


+9.0(1.6) 


+ 10.8(9) 






% 


+9.2(1.3) 


+7.3(1.9) 


+8.0(1.1) 




7T-N 


I = 1/2 






+3.04(66) 

[+2.59] 


[+2.86] 




7 = 3/2 






-1.10(20) 

[-1.29] 


-1.31(22) 

[-1.43] 


7T-7T 


1 = 
1 = 2 






+3.02(17) 

[+3.32] 
-0.924(40) 

[-0.947] 





Figure Captions 



Fig. 1 : R(t) for the nucleon four-point function for the spin singlet and triplet channels 
at f3 — 5.7 and K = 0.160 on an 20 3 x 20 lattice in quenched QCD. Solid line is the linear 
fit to data over 4 < t < 9. 

Fig. 2 : R(t) for the (a) n-N and (b) n-n four-point function at /3 — 5.7 and K = 0.164 
on an 12 3 x 20 lattice in quenched QCD. Solid line is the linear fit to data over 4 < t < 9. 

Fig. 3 : Scattering lengths in lattice units at (3 = 5.7 in quenched QCD as a function of 

m^/rrip. 



1.4 



1.2 - 



1.0 - 



0.8 








8 



12 16 20 



t 



Fig.l 



1.6 



1.4 



1.2 



1.0 



t — i — i — I — i — i — i — I — i — r 



K-N 



1=1/2 

R (t)+0.1 



▲ ▲ 




0.8 - 



0.6 




1=3/2 

h R (t)-0. 1 







1 







8 12 

t 



16 20 



Fig.2(a) 



1.6 



t — i — i — I — i — i — i — I — i — i — r 



t — i — r 



1.4 



1.2 



7I-7C 



R I=0 (t)+0. 1 



▲ ▲ 




1.0 



0.8 - 




0.6 



R I=2 (t)-0.1 



^r» ♦ ♦ ♦ 



1 I I I I I I I I I I I I I I I I I L 







8 



12 16 20 



t 



Fig.2(b) 
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